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Abstract — In this paper, we investigate the problem of non- 
myopic (multi-step ahead) quantizer design for target tracking 
using a wireless sensor network. Adopting the alternative condi- 
tional posterior Cramer-Rao lower bound (A-CPCRLB) as the 
optimization metric, we theoretically show that this problem can 
be temporally decomposed over a certain time window. Based on 
sequential Monte-Carlo methods for tracking, i.e., particle filters, 
we design the local quantizer adaptively by solving a particle- 
based non-linear optimization problem which is well suited for 
the use of interior-point algorithm and easily embedded in the 
filtering process. Simulation results are provided to illustrate the 
effectiveness of our proposed approach. 

I. Introduction 

Wireless sensor networks (WSNs), consisting of a large 
number of spatially distributed sensors, have been used in 
a wide range of promising applications such as battlefield 
surveillance, environment and health monitoring. However, 
due to the limited communication and energy resources, it 
is desirable that only quantized data be transmitted from local 
sensors to the fusion center (FC). 

Quantizer design for target tracking has been recently stud- 
ied in the literature (T]-(4)- In (TJ, the authors employed a 
static quantizer, first proposed in [5], to track a moving target, 
where the optimal quantization thresholds are determined by 
maximizing the Fisher information about the signal amplitude 
contained in quantized data. Although this approach is robust 
and requires minimum prior information about the system, it 
doesn't yield the optimal solution for tracking scenarios where 
the target state is random and dynamic. An adaptive binary 
quantizer and uniform quantizer are proposed in Q and J3J, 
respectively, where the local quantizers are designed at FC 
for every time step. In [4], a more general framework for de- 
signing adaptive identical/non-identical quantizers is presented 
where the trace of the direct conditional posterior Cramer- 
Rao lower bound (D-CPCRLB (6|) is minimized. However, D- 
CPCRLB yields an intractable objective function, which leads 
to a high computational complexity in optimization. Therefore, 
to satisfy the real-time operational requirement of the adaptive 
system, it is essential to seek a tractable objective function and 
develop an efficient algorithm for quantizer design. 

Fisher information matrix (FIM) has been used in (TJ, J3]- 
J5] as the performance metric, where the mean square error 
(MSE) is lower bounded by the inverse of FIM. However, 
the FIM is generally a matrix. It is important to use a suitable 
scalar norm of it to obtain a quantity related to the information 



content. Authors in [4| employ the trace of the inverse of FIM. 
The determinant of FIM is used in (TJ, which is inversely 
proportional to the volume of the uncertainty ellipsoid. In this 
work, we adopt the trace of FIM as the performance criterion 
and will theoretically show that maximizing the trace of FIM 
does not lose the optimality of maximizing FIM in the sense 
of positive semidefinite cone (7), JSJ. Since the trace operator 
is linear, we can also show that the trace of FIM yields a 
tractable objective function for optimization. 

For adaptive quantizers proposed in (2)-(4j, we note that FC 
is required to feed back the quantization thresholds to local 
sensors at every time step (a.k.a., the myopic/greedy design 
strategy J9j). However, continual transmissions might result 
in data collisions and channel congestion. In order to reduce 
the amount of communication, we design adaptive quantizers 
using the non-myopic (i.e., multi-steps ahead) strategy, which 
has drawn recent attention in resource management, e.g., (TJ, 
1 9 1, 1 10). In general, myopic design has lower computational 
complexity than the non-myopic case |9). However, we will 
show that the non-myopic quantizer design can be temporally 
decomposed based on the alternative conditional posterior 
Cramer-Rao lower bound (A-CPCRLB |6|), which, unlike D- 
CPCRLB used in |4], yields a recursive form of the informa- 
tion matrix. With the aid of particle filtering methods |12|, 
the problem of non-myopic quantizer design is expressed as 
a non-linear optimization problem which is easily solved by 
the interior-point algorithm. 

II. Problem Formulation 

In this paper, the task of the WSN is to monitor a sin- 
gle target moving in a two-dimensional Cartesian coordinate 
plane. At sampling time t, the target state is defined by a 
4x1 dimensional vector x t = [x t ,y t , x t ,y t ] where (x u y t ) 
and (x t , y t ) denote the target location and velocity in the 2D 
plane, respectively. The target state evolves according to 



Fx t +w ( , 



(1) 



where w f ~ A/"(0, Q), the state transition matrix F and the 
process noise covariance Q are given by (TJ 
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In A and q denote the sampling interval between adja- 
cent sensor measurements and the process noise parameter, 
respectively. 

We further consider N sensors deployed in a region of 
interest (ROI) and each of them reports a noisy measurement 
in the form of signal power JT], Q 



yl = hl(x t )+vl, ht(x t ) = ^ 



Po 

K) 2 



(3) 



for i = 1,2, ••• , N, where v\ — AT(0,<t 2 ), Po denotes the 
signal power of the source, d\ is the distance between the target 
and the ith sensor, d\ — J (xj — xt) 2 + {yi — yt) 2 , where 
(x^yi) is the position of the ith sensor in the 2D plane. 

Each sensor quantizes its measurement to M bits as given 
below 



-oo < y\ < Yt.i 

L - 1 7 * < y\ < +oo 



(4) 



where u\ denotes the quantized measurement of the zth 
sensor at time step t, L = 2 M , and the vector -y| := 
[yli,"- ,7tL-i\ T corresponds to the quantization strategy 
of sensor i. For notational consistency, let jl = — oo and 
7j L = oo. It is clear from that the probability of a 
particular quantization output I is 



p(ut = J|x t ) = Q( 



(5) 



where Q(-) is the complementary distribution function of 
the standard Gaussian distribution. Under the assumption of 
conditionally independent observations at local sensors, the 
observation likelihood function at time t can be written as 



p(u t |x t ) = IJp(uJ|x t ), 



(6) 



where u t 



,NlT 



denotes the collection of 



quantized measurements from TV sensors. 

A. Alternative conditional posterior Cramer-Rao lower bound 
(A-CPCRLB) 

A conditional posterior Cramer-Rao lower bound (C- 
PCRLB) is proposed in fTT[ by incorporating the history of 
actual sensor observations, which can provide a tighter error 
bound than the conventional PCRLB. Nevertheless, obtaining 
C-PCRLB is not computationally efficient due to the presence 



of the auxiliary Fisher information matrix 1 1 1 Thm. 1]. There 



fore, Zheng et al. in |6| presented an alternative conditional 
PCRLB (A-CPCRLB), which is direct and more compact. In 
this work, we adopt A-CPCRLB as the performance criterion 
for quantizer design. 

Let xo:t and Ui-.t denote the state vector and measurements 
up to time t. Then the conditional mean squared error of 
the state vector xo : t is lower bounded by the inverse of the 
conditional Fisher information matrix (C-FIM) as in (TT1 



Let J(x t+1 |u 1:t ) be the matrix whose inverse equals the 
lower-right corner submatrix of J _1 (xo : t+i|ui : t). Then the 
matrix J(x t+ i|ui : t) provides a lower bound on the mean 
square error (MSE) of estimating x i+ i. As shown in (6] 
Corollary 1], for the linear Gaussian model ([TJ, the C-FIM 
J(x 4+1 |u 1:t ) can be computed as follows, 



(Q + FJ^F 1 



£ P ? +1 {-V££lnp(u t+ i|x t+1 )}, (7) 



where for notional simplicity we use ,] t +i instead of 
J(x t+ i|ui :t ), Vx is the second-order partial derivative with 
respect to x and p1 +1 = p(x t+ i, Ut + i|ui :t ). Note that the first 
term in (j7]i is the prediction of J t using the state evolution 
model and the second term indicates the information based on 
the updated measurements u. t +i at time t + 1. 

B. Non-myopic quantizer design 

For the non-myopic quantizer design, we seek optimal quan- 
tizers defined in Q over the next T w time steps, t+1 : t + T w , 
at time instant t. It is clear from |6]) and |7| that the value of 
C-FIM Jt+ n relies on quantization thresholds of local sensors 
at time t + J] (denoted by -y t+v = hi^,,-" >7t+, ; ]X the 
previous C-FIM Jt+,,-1, and the conditional distribution Pt +7J , 
i.e., p(x t+ ^,u t+7) |ui :t+?) _i), where r\ G {1,2, ••• ,T W } and 
J t is used as prior information. However, during the design 
window t+1 : t + T w , the conditional PDF p c i +v cannot be 
obtained exactly for -q > 1 since the quantized measurements 
u *:i+?;-i are n °t available at time t. Therefore, as in |10) , 
the conditional PDF pl +rj is approximated by its prediction 
p(xj_|_„, Ut_|_Jui : t), which is easily obtained using a particle 
filter (see more details in Sec. |TTT] ). 

To determine optimal thresholds {7t+ ?? } ??= i ! ... t w f° r the 
next T w time steps, we pose the optimization problem as given 
below where we maximize C-FIM at time t + T w , 



maximize J t+Ttu (74+1 , • 



, Jt+T a ) 



subject to Yt+ V ,i < 



, T,„ and i = 



(8) 



.N 



where for notational simplicity we use {■} instead of 
{■} T)= i l ... ,t„> lt+rj is a (L— 1) x N quantizer threshold matrix 
whose element 7| + ; represents the Ith threshold of sensor i 
at time t + rj, L indicates the number of quantization levels, N 
is the number of sensors and T w is the length of time window. 

Note that problem ([8]l is a matrix optimization problem 
which is defined in the positive semidefinite cone [7]. Namely, 
if {it+r)} i s an optimal solution, then for an arbitrary feasible 
solution {jt+ v }, Jt+Tw{{lt +71 }) t Jt+Tw({lt+r,}f\ i.e., the 
matrix Jt+Tw({lt +V }) - J t +Tw({lt+ V }) is positive semidef- 
inite. 

Furthermore, the following Proposition shows that the prob- 
lem ([8]) can be equivalently transformed to T w sub-problems. 
Each of the subproblems has a scalar objective function in 
terms of the trace of the Fisher information matrix with respect 
to the updated measurements. 



E{[XQ:t+l - X():i + l][x :t + l - X 0: t+l] T | Ul : j } > J 1 (x 0: t+l |Ul :t ) . 



A positive/negative semidefinite matrix A is denoted by A y or A ^ 0. 



(9) 



Proposition 1: If problem (|8]l has an optimal solution, then 
the solution of |8]l can be equivalently transformed to the 
solution of T w subproblems, i.e., 

maximize tr \E Pt {-V^lnp(u t+ ,|x t+ ,)} ) 
subject to 7t 8 +I)1 < • • • < 7t+„,z,-i, i = 1) ■ ' ' , N 
for rj = 1,2, ■■ • ,T W , where tr(-) denotes the trace operator 

and Pt+r, ~ P(Xt+r„ U t+7) |Ui :t ). 

Proof: See appendix. ■ 

III. Particle-based non-myopic quantizer design 

In this section, we will show that the problem Q can be 
further expressed in a closed form and solved efficiently with 
the aid of a particle filtering method. 

At time step t + rj, by substituting ([6]) into (|9|, the objective 
function in |9]l can be written as 

N 

E tr (*k* t+ „,„; +0 iu 1:t) {-v*::>pK + ,,ix t+ ,,)}) , 

which indicates that seeking the optimal quantizers of N 
sensors at time t + rj can be obtained by equivalently solving 
a sequence of sub-problems, i.e., 

iphl+r,) = tr (e {-v££M*4J**+*)}) 

7t+jj,l < ' ■ ■ < lt+ri,L-li 



maximize 

Tt+ n 

subject to 



for i = 1,2,- • • ,7V. 

Using the fact that u\ +Tj , x t+)) and ui : < form a Markov 
chain and the identity for standard Fisher information matrix 

CCD 



(ii) 



where x t r and x t j denote the rth entry and jth entry of 
state vector x t , the objective function i/>(7j +r; ) in (lOi can be 
written as 



' d\np{u\\yL t ) 91iyj(i4|x t )' 


= -E 


"9 2 lnp(?tj|x t )' 
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J X t+T|) 



(12) 

where the likelihood p(uJ +t? |x i+?? ) is given by (pi. 

We employ a particle based method to computethe posterior 
PDF. In a SIR filter (H), the posterior PDF p(x t |u 1:t ) is 
approximated by a set of particles {xf ; s = 1, . . . , N s } with 
equal weights 1/N S after the re-sampling process, where N s 
is the total number of particles. Thus, the predicted PDF 



p(x t+)) |ui :t ) in ( 12 1 can be obtained by propagating particles 
xf after rj steps using the state model (JTJ. Then, 



p(x t+r ,\u 1:t ) 



N. 



1 N ' 

=i 



x i+i) 



(13) 



Substituting |5]) and ( p"3j ) into ( [12) , the optimization problem 
TOll can be written as 



L-l 

maximize i]j{~fi +r ,) = E /frt+^Tt+^i+i) 
subject to Y t+7hl < Y t+V ,2 < ■ ■ < 7t+„,i-i 



(14) 



where 

>7t+,,,;+iJ 

and 



ff(7t+?7,ii 7t+7),!+lJ x t+77i r ) 



with 6^ = /i^(xf +J? ) 



L )-Q( 



Tt'+t|,i+l 6 t+ f| 



(15) 



It is clear from (15i that the objective function of (14i is 
non-linear but differentiable. Therefore, the interior-point algo- 
rithm J7) is a well-suited optimization tool for solving problem 



( 14 1. The procedure for non-myopic quantizer design under the 
SIR filtering framework is summarized in Algorithm 1. 



Remark 1: It can be seen from (14i that the complexity 
of the non-myopic quantizer design depends on the number 
of Monte-Carlo particles (i.e., N s ). Fewer number of particles 
would reduce the computation cost but result in worse estima- 
tion performance due to the low accuracy of approximating the 



predicted PDF in ( 13 1. Therefore, it is important to investigate 
the tradeoff between the number of particles and the estimation 
performance. 



(10) Algorithm 1 Adaptive non-myopic quantizer design 



1: At time t, begin with the updated particles xf and weights w* = 

for rj = 1, . . . ,T W do 

Propagate particles by x^ +1 = Fx* + w t 

p(x n+1 \u 1:n _) = Y^s=! 5(x„+i - Xn+l) 



Obtain optimal thresholds for N sensors by solving 

|14) for i = 1, • • • , JV. 
6: end for 

7: Feed {7t+^} 7; =i,... 1 T ro back to local sensors and update par- 
ticles by using the corresponding quantized measurement at 
t+l,...,t + T w . 



A. Binary Quantizers 

For a binary quantizer, i.e., L = 2, the optimization problem 
( |14[ ) becomes unconstrained, i.e., 

maximize /(-co, 7t +J) ,i) + fht+n.u 00 ) , (16) 

Tt+1,1 

whose optimality condition is presented by the following 
proposition. 

Proposition 2: The optimality condition for the binary 
quantizer of sensor i at time t + rj can be expressed by a 
nonlinear equation 



EE 



Oh 



+ ct -EE 



0~X-t+ri,r 

dh\,Jy 



J Q( 7 '+"' 1 



Tt + T),l b t'+r, ■ 



- 2 { 3( jw h; +" )(2Q( 7t ' + ' M h '+ 



+ ■7,1 "1 + .7 ' 



= 0, 



where 7j +J) l represents the quantization threshold of the ith 
sensor at time t + rj, b l ^ v = /iJ +?) (xf +?? ), and h\ +71 (-) is the 
measurement model. 

Proof: The result can be easily obtained by taking the first- 
order derivative of (16) . ■ 

B. Identical Quantizers 



to each other since the target trajectory is almost deterministic 
(and thus predictable) as q — 2.5 x 10~ 3 . 



It is clear from ( 14 1 that local quantizer design relies on the 
sensor location and the predicted measurement h\ +rj {x. s t+rt ), 
which theoretically verifies the statement in [4| that the use of 
identical quantizers at all the sensors leads to performance 
degradation. On the other hand, if we consider a linear 
measurement model y\ = hjx t + v\ for sensor i at time t, 
where hj = ■ ■ • = (e.g., the mean estimation problem 
in fl4)), it can be shown that the problem of quantizer design 
yields identical optimal thresholds for N sensors at every time 
instant due to the effect of identical sensor observation model 
on ([5} and (12). 

IV. Simulation Results 

In our simulations, we consider that N = 9 sensors are 
grid deployed in a 20 x 20 m 2 surveillance area. For target 
motion in we select the sampling interval A = 0.5 
seconds and process noise parameter q £ {0.1,2.5 x 10~ 3 }, 
where the magnitude of q indicates the relative uncertainty 
regarding the target trajectory (see (Tl Fig 3] for an exam- 
ple). The initial state distribution of the target is assumed 
to be Gaussian with mean /jq = [—8.8,-8.8,1.8,1.8] and 
covariance So = diag[cro, er 2 , 0.01, 0.01] where 3<7o = 2. We 
perform target tracking over 10 seconds, i.e., 20 time steps 
(T^ < 20). Sensor measurements are obtained from Q, where 
Po = 1000 and sensor observation noise with a v = 0.1. 
We assume that the fusion center has perfect information 
about the target dynamical model and the noise statistics. 
Observing from simulation results, which are omitted here 
for brevity, using 1000 and 50 particles in target estimation 
and quantizer design, respectively, provides a suitable tradeoff 
between the number of particles and tracking performance, 
which is evaluated in terms of mean square error (MSE) over 
100 trials. 

In Fig.[T[ we demonstrate the tracking performance of the 
2-bits non-myopic quantizer for different time window sizes, 
i.e., T w £ {1,5,10,20}, where the non-myopic design with 
T w = 1 is equivalent to the myopic design, and the non- 
myopic quantizer with T w = 20 becomes an offline quantizer 
since the corresponding A-CPCRLB is calculated offline. For 
comparison, we also present the tracking performance when 
using analog data (AD) and quantized data based on the 
offline Fisher information heuristic quantizer (FIHQ) G). As 
we can see, our proposed quantization strategy yields better 
performance than FIHQ. Specifically, Fig.[T]-(a) shows that the 
estimation performance improves as T w decreases. This is 
because for q = 0.1, the target trajectory has relatively large 
uncertainty so that the accuracy of estimation benefits from 
quantizer design using more sensor measurements. However, 
Fig.[T[(b) shows that the MSE for all values of T w lies close 
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(b) 

Fig. 1. Tracking performance of the non-myopic quantizer with different 
sizes of time window (a) q = 0.1 (b) q = 2.5 X 10~ 3 

In Fig.[2] we present the MSE of temporally identical 
quantizer (I-Quantizer), which refers to the design of quan- 
tizers only for the next time step and then using the same 
quantizers over the entire time window. For comparison, the 
MSE of nonidentical quantizer (N-Quantizer) with T w — 20 
is also plotted. Simulation results show that the I-Quantizer 
yields worse performance than the N-Quantizer even with 
a small window size (i.e., T w = 2). This is because in a 
tracking scenario the target state is random and dynamic, 
which leads to a large innovation error by using temporally- 
identical quantizers, although the identical design can save 
energy and computation cost. 

V. Conclusion 

In this paper, we considered the problem of target tracking 
with quantized data in a WSN, where the optimal local quan- 
tizers are determined using a non-myopic strategy. Using the 
alternative conditional posterior Cramer-Rao lower bound (A- 
CPCRLB) as the performance metric, we theoretically showed 
that the non-myopic quantizer can be designed separately 
for each time instant. With the help of a particle filtering 
method, this problem can be expressed in a closed form 
and solved via the interior-point algorithm. Simulation results 
demonstrated the effectiveness of our proposed approach. In 
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Fig. 2. MSE performances of 2-bits quantizer as q = 2.5 X 10 



the future, we will consider the effect of channel statistics on 
quantizer design. We will also consider a unified non-myopic 
optimization framework for resource management problems 
such as sensor selection and bit allocation. 

Appendix 

Proof of Proposition 1: According to (T7}, problem ([H} can 
be decomposed into two subproblems 



maximize 



[-V^ : in P (u t+T Jx i+T J} 



and 



subject to 7j +T i < • • • < f t+Tn ,L- 



maximize (Q + FJ t } T _ 1 F T ) 
{~tt+ q } 

subject to T | +I)il < ■ ■ ■ < 7 t +„,z,-i 
77 = 1, • - • , T„ — 1 and i 



(17) 



,N 



where Q and F are given by the process model ([T]). 

Note that Q + FJ^r T _ 1 F T is positive definite since Q 
is positive definite and the information matrix Jt+T m -i is 
positive definite (where we assume it is invertible). Then, 
problem ( fTT] ) can be written as 



minimize Q + FJ t 

subject to ji +rl>1 < 
?7 = 1 , ■ ■ 



t+T„-l 



F T 



(18) 



,N 



where we use the fact that, for any positive definite matrix, if 
A y B then B^ 1 > A~ l . 



Since F is invertible, the problem ( 18 1 is equivalent to 



maximize J t+ t. iv - 1 

subject to yl +rtA <■■■< 7l +v ,L-i 

tj =!,••• , T w — 1 and i = 1, 



,N 



Similarly, after T w recursive decompositions, the problem 
dH) can be decomposed into T w sub-problems given by 



E Pt+n {-v£:>p( Ut ^ix t+ „)} 



subject to 7| +J?1 < • • • < 7t +I) ,L-i, i = 1, 



maximize 

74+r, 



,iV 



(19) 



Then by |8j Lemma 3.1], the problem ( 19 1 is equivalent to 
the problem (|9|, where for clarity , we reiterate the |8] Lemma 
3.1] as below. 



Consider two optimization problems 



M(x) 



max 
xeS 



max tr(A/(x)) 

xG5 



G4i) 
(A 2 ) 



where M(x) is a matrix for an arbitrary x € S, S specifies the 
constraint on x. If the problem {A~i\ has an optimal solution, 



then the problem §Ai\ is equivalent to CA 



The proof of |8] Lemma 3.1] includes two parts. First it can 
be shown that if xi is the optimal solution of \A\) , then for 
arbitrary x G S, M(x) < M(xi) which yields tr(Af(x)) < 
tr(M (xi)). Thus, Xi is also the optimal solution of (A% 1. On 
the other hand, if x 2 is the optimal solution of (A~ 2 



then we 

have tr(M(xi)) < tr(M(x2)). Note xi is the optimal solution 
of §A^ which implies tr(M(x 2 )) < tr(M(xi)). Thus, we 
obtain tr(M(xi) - M(x 2 )) = 0. By tr(M(xi) - M(x 2 )) = 
and M(xi) - Af(x 2 ) y 0, we have M(xi) = M(x 2 ). 
Therefore, x 2 is also the optimal solution of \A X \ . ■ 
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